Effect of electron-electron interaction on the Fermi surface topology of doped 

graphene. 
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The electron-electron interactions effects on the shape of the Fermi surface of doped graphene are 
investigated. The actual discrete nature of the lattice is fully taken into account. A 7r-band tight- 
binding model, with nearest-neighbor hopping integrals, is considered. We calculate the self-energy 
corrections at zero temperature. Long and short range Coulomb interactions are included. The 
exchange seff-energy corrections for graphene preserve the trigonal warping of the Fermi surface 
topology, although rounding the triangular shape. The band velocity is renormalized to higher 
value. Corrections induced by a local Coulomb interaction, calculated by second order perturbation 
theory, do deform anisotropically the Fermi surface shape. Results are compared to experimental 
observations and to other theoretical results. 
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I. INTRODUCTION 

Since its discover}^^!^ graphene, a two-dimensional 
(2D) single crystal thermodynamically stable formed by 
a single-layer of carbon atoms ordered in a honeycomb 
lattice, has been thoroughly investigated. It forms the 
basic block of carbon nanotubes, fuUerenes, graphite, 
and graphite intercalation compounds. The 2D electronic 
properties are well described by a 7r-band tight-binding 
(TB) model^. The valence and conduction 7r-bands touch 
only at the six corners, K, of the 2D Brillouin zone 
(BZ); the degeneracy point of the valence and conduc- 
tion bands is often termed Dirac point. At half-filling, 
undoped graphene, the Fermi level Ep lies at the Dirac 
point. The low-energy physics of a perfect graphene sheet 
is described by the relativistic Dirac equation. The dis- 
persion relation turns up to be isotropic and linear near 
Ep. The low-energy excitations of the system are Dirac 
fermions with zero effective mass and a vanishing density 
of states at the K points. Because of these peculiarities, 
graphene is considered a model system to investigate ba- 
sic questions of quantum mechanics. Due to its transport 
properties graphene is a promising material for nanoelec- 
tronic applications^. 

Improvements in experimental resolution have led to 
high precision measurements of the Fermi surface (FS), 
and also to the extraction of the many-body effects from 
the spectral function, as reported by angle resolved pho- 
toemission spectroscopy (ARPES) experiments^ in dif- 
ferent materials. The recent isolation of an atomic layer 
of graphite, graphene, has renewed the interest in the 
physics of the three-dimensional (3D) graphite and new 
aspects of the electronic properties of 3D graphite have 
been observed with improved experimental techniques^. 

In a high resolution ARPES study of disordered 
graphite samples, coexistence of sharp quasiparticle dis- 
persions and disordered features was found^, and was ex- 
plained in terms of Van Hove singularities (VHS) in the 



angular density of states. Later on, by using ARPES, 
the linear and isotropic dispersion of the bands, near the 
three-dimensional BZ corners (H points) of graphite, has 
been directly observed coexisting with parabolic disper- 
sion bands^. The constant energy maps taken near the 
H point present circular shape from Ep to — 0.6cV. This 
circular shape combined with the linear dispersion found 
near the BZ corners H suggests that the dispersion shows 
a cone- like behavior near each point H, similar to that 
expected for graphene. The constant energy maps start 
to deviate from the circular shape at — 0.9eV, and at 
— 1.2eV a rounded triangular shape is observed. A lin- 
ear energy dependence of the quasiparticle life-time has 
been as well measured by ultrahigh resolution ARPES on 
high quality crystals of graphitei^. The low-energy exci- 
tations seem to be dominated by phonons, while those for 
higher energies arc characterized by the electron-hole pair 
creationiii. The quasiparticle life-time had been studied 
as well by ultrafast time resolved photoemission spec- 
troscopy (TRPES}ii and a linear uj dependence was re- 
ported. Anisotropy of quasiparticle life-times has been 
reported by TRPES in highly oriented pyrolytic graphite 
(HOPG) as well as an anomaly in the energy dependence 
between 1.1 and 1.5eV, in the vicinity of a saddle point 
in the graphite band structure^^. 

A linear energy dependence of the quasiparticle life- 
time had been theoretically predicted for graphiteiS. ne- 
glecting the interlayer hopping. Even including an in- 
tcrlayer hoping of 0.25eV, an anomalous quasiparticle 
life-time was obtained in graphite^ within the GW ap- 
proximation, with a linear energy dependence along the 
Kr direction for energies well above the interlayer hop- 
ing. A discussion about the interlayer coupling strength 
in graphite can be found in RefSfii^^, along side with 
an overview of recent experiments in both graphene and 
graphite samples. 

In a combined ARPES and theoretical ab-initio quasi- 
particle study of the 7r-band structure and the Fermi sur- 
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face in graphite single crystals, it is found that electron- 
electron correlation plays an important role in semi- 
metallic graphite and should be taken into account for the 
interpretation of experimental results^^. The electronic 
correlations renormalize the electronic dispersion increas- 
ing the Fermi velocity. The equi-energy contours of the 
photocmission intensity show trigonal warping (higher by 
a factor of about 1.5 if compared to graphcnc) around 
the KH direction of the graphite 3D BZ, in both Lo- 
cal Density Approximation (LDA) and TB-GW calcula- 
tions, from — O.leV. Correlation effects are found to be 
stronger as the energy increases and differences between 
LDA and GW results are more noticeable at — 0.4eV than 
at -O.leVi^. 

The recently available 2D graphene samples have been 
intensively investigated. From the experimental point of 
view graphene presents advantages with respect to other 
2D systems. Graphene can be controlled externally and 
exposed to vacuum therefore can be directly probed by 
different techniquesii^i^. Electron-electron interactions 
in graphene are expected to play an important role due 
to its low dimensionality and many-body effects have re- 
ceived great attention. The quasiparticle dynamics in 
graphene samples has been addressed by high-resolution 
ARPESii. It was found that the conical bands are dis- 
torted due to many-body interactions, which renormal- 
ize the band velocity and the Dirac crossing energy Ejj. 
Electron-hole pair generation effects are important near 
the Fermi energy and electron-phonon coupling contribu- 
tion to the self-energy is also important in the Fermi level 
region: an electron-phonon coupling constant of A « 0.3 
is deduced with the standard formalismii. Around E]j, 
electron-plasmon coupling is invoqued to explain the 
peak in the imaginary part of the self-energy ImS, found 
just below Ed, whose width (and intensity) scales with 
E]j. Although the three scattering mechanisms con- 
tribute to strongly renormalize the bands, in Refiii it is 
claimed that the quasiparticle picture is valid over a spec- 
tacularly wide energy range in graphene. More recently, 
the doping dependence of graphene electronic structure 
has been investigated by ARPES in graphene samples 
at different dopingsi^. Upon doping with electrons, the 
Fermi surface grows in size and deviates from the circular 
shape showing the trigonal warping, finally evolving into 
a concave triangular shape. An electron-phonon coupling 
to in-plane optical vibrations is proposed to explain the 
experimental results. The electron-phonon coupling con- 
stant, extracted from the data, presents a strong depen- 
dence on k, with maximal value along the KM direction. 
The presence of a Van Hove singularity in the KMK 
direction, confirmed upon doping, could be a possible 
explanation of the enhancement of the electron-phonon 
couplingi^. 

The layered nature of graphite/graphene as well as the 
presence of the VHS in the density of states near the Ep, 
which would enhance many-body effects, make contact 
with the physics of cuprates superconductors. The sim- 
ilarities between graphene and the cuprates have been 



already noticedi^ and the important role of many-body 
effects in the basic physics of graphene has been investi- 
gated earlier—. 

ARPES investigation of graphene samples epitaxially 
grown on SiC substrate, have reported the observation of 
an energy gap of « 0.26eV at the K point^i. It is pro- 
posed that the opening of the gap is induced by the in- 
teraction with the substrate on which graphene is grown, 
that breaks the A and B sublattice symmetry^i. 

From the theoretical point of view, graphene offers 
many possibilities and the electron-electron interaction 
induced effects have been widely analyzed. The self- 
energy have been object of special interest, since it gives 
relevant information from fundamental properties. Fur- 
thermore, theoretical results can be compared to re- 
cent ARPES reported self-energies. Many approxima- 
tions have been followed to investigate the self-energy in 
graphene. The inelastic quasiparticle lifetimes have been 
obtained within the GqW approximation^ with a full dy- 
namically screened Coulomb interaction. The scattering 
rates calculated for different carrier concentrations are in 
good agreement with ARPES data from Bostwick et alJ^ 
without including phonon effects, contrary to the exper- 
imental interpretation of the dataF^. The nature of the 
undoped and doped graphene has been as well discussed 
theoretically in terms of the behavior of the imaginary 
part of the self-energy^^: a Fermi liquid behavior is found 
for doped graphene while the zero doping case exhibits 
a quasiparticle lifetime linear and a zero renormalization 
factor indicating that, close to the Dirac point, undoped 
graphene behaves as marginal Fermi liquid, in agreement 
with previous theoretical worki^^^l. 

By evaluating exchange and random-phase- 
approximation (RPA) correlation energies, an en- 
hancement of the quasiparticle velocities near the 
Dirac point is found in lightly doped graphene taking 
into account the eigenstate chirality^^. The role of 
electron-electron interactions in the ARPES spectra 
of a n-doped graphene sheet has been theoretically 
investigated^^ by evaluating the self-energy within the 
RPA and turned out to be important when interpreting 
experimental data. Recently^i^ the validity of the 
RPA in the calculation of graphene self-energy has been 
discussed. The approximation is shown to be valid and 
controlled for doped graphene where the Fermi level is 
shifted up or down from the Dirac point. The RPA fails 
for undoped graphene, where Ep lies at E]j. In most 
of the self-energy studies graphene is described by the 
massless Dirac equation in the continuum limit. 

In this work, we calculate the corrections induced by 
the electron-electron interaction on the electronic band 
structure. We focus on the corrections to the Fermi sur- 
face shape of doped graphene. The Fermi surface is one 
of the key features needed to understand the physical 
properties of a material and its shape provides impor- 
tant information. Due to the 2D character of graphene, 
electronic interaction should be important. We consider 
the TT-band tight-binding model taking into account the 
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discrete nature of the honeycomb lattice in order to in- 
vestigate doped graphene and possible correlation effects 
in the trigonal warped topology of the Fermi surface. We 
calculate first the exchange self-energy considering long- 
and short-range Coulomb interactions. Besides the renor- 
malization of the band velocities, a deformation of the 
trigonal warped Fermi surface is found at this level. The 
second-order self-energy induced by an onsite Coulomb 
interaction is as well calculated. The corrections to the 
Fermi surface shape are found to be anisotropic. 

The paper is organized as follows. In section II the 
model is presented and we explain the self-energy cal- 
culation method. Section III presents the results of the 
calculation and Section IV contains a discussion of the 
results compared to experimental data and to other the- 
oretical results and the main conclusions of the work. 



II. THE METHOD 
A. The model for the graphene layer 



FIG. 2: Constant energy contours obtained with Eq.Q. The 
thick black lines corresponds to the Van Hove filling. 




Graphene is an atomic layer of carbon atoms arranged 
in a honeycomb lattice with two atoms per unit cell, as 
shown in Fig 1(a) The distance between nearest neigh- 



bor atoms is a = 1.42 A, and the primitive lattice vectors 
are Si and a.2- The Brillouin zone is an hexagon, as de- 
picted in Fig 1(b) We adopt the 7r-band tight-binding 



model with only nearest-neighbor hopping^, since it cap- 
tures the main physics of the system as probed by more 
realistic models and by experimental results. 
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taken as zero. The operator a\ ^ i^i.cr) is the creation 

(annihilation) operator of an electron at site on sub- 
lattice A with spin a (cr =1, |) (an equivalent definition 
corresponds to sublatticc B). By Fourier transformation 
to the momentum space we have 



Ti-kin (k) — t 



5(k) 
S*(k) 



(2) 



The function .g(k) 



■^/3, 



g{k) = - e-'^y'' + 2e^*=«"/2 cos I ^k,.,a 



(3) 



is the structure factor of the honeycomb lattice. Diag- 
tializing Eq.Q the 
A<:|(?(k)| is obtained 



onalizing Eq.Q the energy dispersion relation e°(k) 



FIG. 1: (a) Schematic representation of the structure of the 
honeycomb lattice, open and solid points represent the two 
inequivalent atoms. The basic vectors Si and §2 of the lattice 
are shown. The three fi connect nearest neighbor atoms, (b) 
First Brillouin zone of the honeycomb lattice. 

The kinetic term of the Hamiltonian, considering only 
nearest-neighbor hopping, will be (we choose h = 1 
throughout this paper) 



=5(£)=At 



l+4cos( -aky IcosI 



'^/3 



,ka: +4C0S2| 



(4) 

where A = —{+) for the valence (conduction) band. 
The two bands are degenerated at the six corners of the 
BZ, K points. The corresponding Bloch wave functions 
are 



-Hk.n = E ^A'^ + ^-^^ (1) ^^.^iv) = i= ( /C^g(r) + A-^/C,£(f) 1 (5) 

(y).'^ ' v2 \ |5(k)| / 

where t = 2.82eV is the nearest-neighbor hopping pa- The tight-binding functions arc built from the atomic 

rameter, and the site energy of the 2pz atomic orbital is 2pz orbitals^S, (/)^(r) 
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(6) 



where A and 5 stand for the area of the unit cell and 
the crystal respectively; p„ = nisii + n2a2 is the lat- 
tice vector, and fi^2 define the position of the two in- 
equivalent atoms in the unit cell. We use Ti = and 



T2 = i(ai -I- (sec Fig 1(a)). At half-filling (undoped 



graphene) the Fermi energy lies at the common point 
of the two bands (we take this energy as our zero en- 
ergy) and the Fermi surface is formed by six points at 
the six BZ corners. The constant energy contours for the 
dispersion relation Eq.(I3]) arc depicted in Figl^l These 
six isolated points (only two of them are inequivalent) 
are known as Dirac points because around them, by a 
long-wavelength expansion, the kinetic energy term of 
the Hamiltonian can be approximated by the 2D Dirac 
equation for massless fermions. Upon doping, by follow- 
ing the constant energy maps shown in Fig|21 the FS 
points develop into circles and, eventually, the FS adopt 
the rounded triangular shapes . As can be seeing in Figj^l 
around each of the six BZ corners the energy lines arc the 
same but rotated with respect each other. Therefore, in 
the following, we will show the results in the vicinity of 
one of the corners of the BZ. 

The electron interaction term of the Hamiltonian Tiint 
includes the Coulomb interaction e^/eo|fy|, where |ry | 
indicate distances between sites of the honeycomb lat- 
tice and £0 is the dielectric constant. We will study the 
two limiting cases, long- and short-range Coulomb repul- 
sion. Among the short-range interactions, we consider 
interactions between electrons on neighbor atoms of the 
honeycomb lattice. 

We also analyze corrections due to an on-site Coulomb 
repulsion U between electrons with opposite spin on the 
same Pz atomic orbital. Although graphene is consid- 
ered a weakly correlated system, the effects of an on-site 
Coulomb interaction, Hubbard term, on the electronic 
properties of the honeycomb lattice have been investi- 
gated in different scenarios. The phase diagram of the 



Hubbard model in the honeycomb lattice has been stud- 
ied by a variety of techniques. Different instabilities 
have been found, from the Mott-Hubbard transition^, 
charge- and spin-density wave^i*^, to superconductivity 
and magnetic phases^^i^i^^. In our calculation we con- 
sider small values of the on site repulsion, of the order of 
the hopping parameter for graphene. 



B. Calculation of the self-energy 

1. Long- and short-range Coulomb interaction: exchange 
self-energy. 

We calculate first the exchange self-energy contribu- 
tion, that corresponds to the one-loop diagram shown in 
FiglHJa). The self-energy has the form: 



S^(k, iujn)= ^ '^^+^' ^^n+ii^n)vxx' (k, k+q) 

A' = ± q w„ 

(7) 

where (3 = l/ksT, ujn and Vn are fermionic and bosonic 
Matsubara frequencies respectively, and G'5(k, iujn) is the 
bare single-particle Green's function of an electron with 
momentum k and band-index A. The Coulomb interac- 
tion matrix elements vxy (k, k + q) between states |k, A) 
and |k -I- q. A'), are given by 




(a) 



FIG. 3: (a) First-order self-energy diagram corresponding to 
the exchange contribution, (b) Second-order self-energy dia- 
gram. 



i;,,v(k,k + q) = / dfi / d?(*k.;,(fi)|(*k+q,A'(?i -?)|-|*k.A(fi -?)>l*k+a.A'(ri)> 



= F(q)FAA'(k,k + q) 



k.AV 



k+q,A' 



(8) 



r 



where V{i\) is the 2D Fourier transform of the Coulomb 
interaction 



F'^(q) 



27re^ 
eo|q| 



(9) 



We are also interested in the effects of short-range 
Coulomb interactions between electrons on neighbor 
atoms of the honeycomb lattice. The Coulomb interac- 
tion in momentum space between neighbor sites (matrix 
elements), described in Appendix 1X1 can be expressed as 
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y^'-(q) = V\2 cos{qya) + 4 cos 



V3 



-qxd COS 



-qya 



(10) 

We take V = t for the nearest-neighbor interaction 
strength. The function FxA'(k, k') in Eq.® arises from 
the overlap of the wave functions obtained by diagonal- 
izing the 7r-band tight-binding Hamiltonia n^^i"^^'^^ , 



FAA'(k,k + q) 



^/^(|q|) 



1 + AA 



,g(k)g*(k + q) 



|g(k)5(k + q) 



(11) 



Therefore the correct symmetry of the lattice is in- 
cluded in UAA'(k, k -I- q). /(q) comes from the matrix 
elements that contain the 2pz wave function of carbon 
atoms (pzif) and can be approximated by the unitj^. 
The function FAA'(k, k + q), close to the Dirac points, 
reduces to the form, 



Faa' (k, k + q) = (1 + AA'cos% £^.)/2 (12) 

where 9^ is the angle between k and k + q. 
The form of the sublattice overlap matrix element for 
graphenc, _F'AA'(k, k + q), given in Eq. p2p appears in 
the theoretical studies based on graphenc masslcss Dirac 
equation continuum mode P^i^^'^^i'^^'^^ . 

Considering the T = limit, after performing the sum- 
mation of Matsubara frequencies, we can write 



dq 



5]-(k)=-^ ^V{q)Fx,yih,k+^M^^-el,ih+ci)) 

(13) 

where & is the Heaviside unit step function. The cal- 
culation of S'J (k) requires a momentum integral over the 
first BZ. Note that, since we are considering the dis- 
creteness of the lattice, the q-integral does not have the 
logarithmic ultraviolet divergences appearing in the con- 
tinuum model. To carry out the integral in momentum 
space we divide the hexagonal BZ in two regions: the 
central region with small momenta q, where the contri- 
bution of the long-range interaction is important, and the 
rest of the BZ with larger momenta, where short-range 
interactions, namely, nearest-neighbor interactions, are 
dominant. In the boundary between both regions the 
potential functions are smoothly matched. Details are 
given in Appendix B. As stated above, the matrix ele- 
ments i^AA' corresponding to intraband (A = A' = + and 
A = A' = — ) and intcrband (A = + and A' = — ) exci- 
tations are calculated by Eg. pTj) . Finally the correction 
of the dispersion relation due to the interactions can be 
calculated from: 



where e°(k) was given in Eq.Q. 

^From Eg. p^ we obtain the corrections to the FS 
shape for finite values of the doping. Notice that we 
address the corrections to the conduction band, taking 
into account the lattice symmetry. 



2. Local Coulomb interaction: second-order self-energy. 

In order to study the effects induced by local in- 
teractions in the Fermi surface topology we compute 
the second-order self-energy. There are two diagrams 
that renormalize the one-particle Green's function up to 
second-order in perturbation theory. The Hartree dia- 
gram gives a contribution independent of momentum and 
energy and, hence, docs not deform the FS. The two-loop 
diagram, depicted in Figl3]^b), modifies the FS topology 
through its k dependence and changes the quasiparticle 
weight through its uj dependence. 

We follow here the method explained in ref4S, study- 
ing the interplay between the electron correlation and the 
FS topology. We assume that the effect of high-energy 
excitations on the quasiparticles near the FS is integrated 
out, leading to a renormalization of the parameters of the 
Hamiltonian. Since we are interested in low-temperature 
and low-energy processes, only the particle-hole excita- 
tions within the energy scale about the Fermi line defined 
by the cutoff A, are taken into account. We compute the 
second order self-energy assuming that a FS, dressed by 
the corrections due to the momentum independent inter- 
action U, can be defined. We consider electron doped 
graphene, thus the FS lies at the conduction band. No 
singularities are reached for the values of the doping con- 
sidered. 

We will perform an analytical calculation of the 
second-order self-energy. We begin by momentum ex- 
panding the function (7(k) given in Eq.([3]), around two 
inequivalent Dirac points, obtaining 



ff(k) 



a{kx + iky) — {kx — iky) around K 



- ikyY around K' 



(15) 

where K = (47r/3V3a,0) and K' (-47r/3V3a, O) . 
After diagonalizing the non-interacting Hamiltonian, we 
find the simplified dispersion relation 



where again A = ± for the conduction and valence 
band, v = ± for valley K and K' respectively, and 
k = |k|. To second order in perturbation theory the 
renormalized FS is given by the solution of the equation: 



6A(k) = e^(k) + SX(k) 



(14) 



e?.(k)-ReSi'l(k,^ = 0) = 



(17) 



6 



where ReEj^ ^(k, is the real part of the self-energy of 
an electron of the A band with momentum k and valley 

f 21 ^ 

index v. The Rel]\j,(k, w) for a momentum indepen- 
dent interaction as the Hubbard U has been computed 
from the imaginary part of the self-energy following the 
method explained in^ and turns up to be, 



ReI]A,^(k,w = 0) 



3 [/^g^ A^sgn(&(k)) 
32^7r3 4(k)&(k) 



(18) 



where the frequency integral has been restricted to the 
interval < a; < A, A being the high energy cutoff taken 
of the order oit. The Fermi velocity vpi^) and the curva- 
ture &(k) of the non interacting FS are given by the first 
and second derivative of ^ (k) with respect to the mo- 
mentum, respectively. Both, u_F(k) and 6(k), are related 
to the parameters of the Hamiltonian. For consistency 
we consider a weak local interaction [/ < A, below the 
energy cutoff. Forward and backward scattering channels 
are considered in the calculation. As explained in^, this 
method of calculation of the self-energy corrections does 
not depend on the microscopic model used to obtain the 
electronic structure and the FS. Simple analytical ex- 
pressions of the effects induced by the interactions are 
deduced from local features of the Fermi surface. 



III. RESULTS 

A. Corrections induced by the exchange self-energy 

We now explain the numerical results obtained for the 
exchange self-energy calculated from the formula Eq.([T3]). 
As stated above, in the momentum integration we con- 
sider the long-range Coulomb interaction in a region of 
the BZ around F, and the nearest-neighbor interaction in 
the rest of the BZ. The results show that the self-energy 
gives the stronger corrections at the boundary lines of 
the BZ, keeping the symmetry of the lattice. The cor- 
rections to the FS topology are found to be small, but 
not negligible, as can be observed in Fig 4(a) for a doping 
density of n ~ 5 x 10^^ electrons per cm corresponding^ 
to a /X = 2.4eV. The self-energy corrections enhance the 



curvature of the sides and round the vertices of the trian- 
gular FS. This result agrees with the renormalization of 
the warping term found in Ref4^ within one-loop renor- 
malization group. One of the conclusions reached in^ 
is that the Coulomb interaction tends to suppress the 
warping term making the energy surfaces more isotropic. 

The self-energy effects are more noticeable in the band 
slope, as can be observed in Fig |4(b)[ where the disper- 
sion relation, calculated without and including the self- 
energy corrections, is shown. The exchange self-energy 
enhances the velocity of the bands by a 20% for the used 
parameter values {V = t — 2.82eV), renormalizing the ki- 
netic energy. The enhancement of the velocity has been 




FIG. 4: Color online, (a) Non-interacting (blue) and interact- 
ing (red) Fermi lines around the K point for a doped graphene 
layer with ^ « 2.4eV due to the exchange self-energy. The 
origin of the coordinate system is at the V point of the BZ. 
(b) Bare (blue) and renormalized (red) bands for ^ = 0.5eV. 



obtained^^i^i^i^^ previously. Calculating the self-energy 
within the on-shell approximation, a renormalization of 
the the velocity compatible with Fermi liquid behavior 
was obtained for doped grapheno2^. The increase of the 
velocity in lightly-doped graphene has been attributed 
to the loss in exchange energy when crossing the Dirac 
point, switching the quasiparticle chirality, by evaluating 
exchange and RPA correlation energie a^^i^^ . All these 
theoretical works are based on the massless Dirac model 
for the low-energy excitations of graphene. We find that 
the exchange energy effects on the FS topology are small 
for doped graphene, in good agreement with the ARPES 
results about the evolution of the graphene FS shape with 
doping^ii^. 



B. Corrections induced by a local interaction 

In this section we analyze the corrections induced by 
a local interaction. The real part of the self-energy is 
computed from Eg. p^ which gives the second order 
perturbation theory renormalization of the Green's func- 
tion. As explained in Ref4ii this electron-electron inter- 
action self-energy depends on local features of the non- 
interacting FS, as the Fermi velocity and the curvature of 
the Fermi line. We limit ourselves to the weak coupling 
regime as seems generally accepted for graphene, and 
where the perturbation approach is justified. In Figl5]the 
Fermi surfaces, at two different doping levels, have been 
represented around the corner K of the BZ. It should be 
notice that the K' counterpart has to be considered, (see 
Figl2] where the hexagonal BZ is represented) in order to 
include all the possible scattering channels. 

At low doping, for /i w 0.2eV corresponding to a 
doping density of about 3.9 x 10^^ electrons per cm^, 
the FS has a circular shape and the correction found is 
small but appreciable. The correction depends on k, as 



shown in Fig 5(a) By increasing the doping level the FS 
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FIG. 5: Color online, (a) Non-interacting (blue) and interact- 
ing (red) Fermi lines around the K point for a doped grapliene 
layer with fj, ~ 0.2eV due to a local Coulomb repulsion, (b) 
Same as (a) but for fi « 2.4 eV. Notice the diflerent scales of 
the coordinate axes. The origin of the coordinate system for 
this case is placed at the K point. 



adopts the round triangular shape. At a doping density of 
71 « 6 X 10^^ per cm^, corresponding to a chemical poten- 
tial of /i « 2.4eV, the effects of the self-energy are notice- 
able, and present a strong anisotropy. The deformation of 
the Fermi surface is maximum along the KM direction 
of the BZ, as can be observed in Fig |5(b)[ The topol- 
ogy of the FS is changed by the correction to a concave 
shape. The underlying hexagonal symmetry (27r/3 rota- 
tion around F point) is preserved. This correction, with 
a maximum along the KM direction, is consistent with 
the Renormalization Group analysis in^, which shows 
that the Van Hove point defines a stable fixed point, so 
that interactions deform the Fermi surface towards the 
saddle point in the band dispersion. Surprising enough, 
this deformation, induced by pure electron-electron in- 
teraction self-energy, presents a strong similarity with 
the deformation found by ARPESi^ in graphene. On 
Fig. 3 (a) from Refi^^ the plot of the Fermi contours de- 
rived from curvefitting the data, obtained at various dop- 
ings, are presented. The anisotropic deformation of the 
Fermi contours are attributed to electron-phonon cou- 
pling to the graphene in-plane optical vibrations. The 
electron-phonon coupling constant Att (k) extracted from 
experimental data shows a strong anisotropy (around 5:1 
ratio at highest doping) and a much higher strength than 
would be expected for 7r-bands and optical phonons in 
graphene. It is argued that the abruptness of the kink 
(change of the band velocity) and the broadening of the 
bands suggest a strong mass renormalization by electron- 
phonon coupling. The divergence of (k) along the KM 
direction is explained as due to the VHS along the KMK 
boundary line of the BZ, which is reached upon further 
dopingi^. The discrepancy between the experimental 
and theoretical electron-phonon coupling constants has 
been addressed in ref.— by using electron-phonon matrix 
elements extracted from density functional theory sim- 
ulations. It was found that, including finite resolution 
effects, the discrepancy between theory and experiment 



was reduced from a factor of 5.5 to 2.2. 

The presence of the strong VHS reveals some similar- 
ity to the cuprates^. The similarities between graphene 
and the cuprates have been already remarked by Bena 
and Kivelsoni^ not only because of the VHS but as 
well because of the similarity between the band struc- 
ture of graphite and that of the nodal quasiparticlcs in 
the cuprate superconductors. The existence of an ex- 
tended VHS, as probed by the flatness of the bands in 
both graphene and cuprates, is important because of its 
effects on the divergences of the density of states. The 
VHS is related with the interactions which can be en- 
hanced in its proximity and could increase the supercon- 
ducting critical temperature^^, in both cuprates and in 
graphite intercalated compounds, when Ep is placed at 
the VHS. 

In our calculated self-energy, considering only a weak 
local electron interaction, the proximity to the VHS en- 
hances its correction and it reaches its maximum at the 
Van Hove band filling, where the Fermi level reaches the 
VHS {fj, = t). Therefore, at low doping, when the FS 
lies very close to the K point, the corrections to the FS 
shape are very small since the circular line lies far from 
the VHS location. On the other hand, at higher doping 
the Fermi line lies close to the VHS, specially the vertices 
of the triangle, and therefore the deformation is higher in 
these directions. In RefJ^ by studying the evolution of 
the coupling parameter A^(k) with the doping level it is 
found that, while the minimum coupling strength grows 
only slowly with doping, the maximum coupling param- 
eter diverges as the corresponding segment of the Fermi 
contour approaches the VHS at the M point. 

The corrections induced by the self-energy calculated 
from Eg. p^ depend on local features of the non- 
interacting FS such as the Fermi velocity vp and the 
curvature of the Fermi line. The main self-energy effects 
in the square lattice were found to occur when the FS 
reaches the border of the BZ, where the real part of the 
self-energy presents pronounced dips close to the saddle 
points because the vp vanishes at the VHS^. The same 
behavior is found here for the honeycomb lattice. 

Because of the relation between the real and imaginary 
part of the self-energy, a large self-energy correction to 
the FS implies a large quasiparticle decay rate in the 
vicinity that region. 



IV. DISCUSSION AND CONCLUDING 
REMARKS 

We have calculated the corrections induced by 
electron-electron interactions in the band structure of 
graphene. We limit ourselves to the weak coupling 
regime, as it is generally assumed to be appropriate 
for graphene. At first order and zero temperature, the 
exchange self-energy induces a deformation which con- 
trarest the trigonal warping of the Fermi surface of 
moderately doped graphene, in agreement with one-loop 
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renormalization group calculations^. The sclf-cncrgy 
corrections round the triangular shape. We find that the 
effect of the long-range Coulomb repulsion is stronger at 
low dopings, when the Fermi level lies near the Dirac 
point. Due to the circular shape of the FS at these fill- 
ings, the corrections to the Fermi surface topology are 
small and only a renormalization of the chemical poten- 
tial is relevant at this level. On the other hand, the Fermi 
velocity is enhanced. 

The self-energy calculated by second-order perturba- 
tion theory, considering a local Hubbard interaction, 
shows a strong k-dependence, therefore the deformation 
induced in the FS topology is anisotropic. The defor- 
mation is highest in the KM direction of the BZ, and 
the maximum value is reached near the Van Hove filling, 
when the Fermi level is very close to the VHS (M point 
in Fig 1 (b) ) . This deformation of the FS shape appears 



to be consistent with the anisotropic deformation found 
by ARPES attributed to a coupling to the graphene in- 
plane optical phononsi^. A different origin is proposed 
in^, where the large and anisotropic values of the ap- 
parent electron-phonon coupling measured by ARPES in 
graphene samples doped with K and Ca— , are explained 
by the strong non-linearity in one of the tt* bands below 
Ei? due to its hybridization with the dopant Ca atoms^. 
In the present work, it is induced by a purely electronic 
interaction. The difficulties in the determination of the 
origin of some features of the ARPES spectra have been 
discussed on a theoretical basis in ref4i. 

We now comment on the connection between our 
work and recent theoretical work on correlation effects 
on graphene. Most of the calculations of the self- 
gj^gj.gjj23j24|25j26 j^^ye been carried out for undoped or 
lightly doped graphene, in the continuum limit, where 
the conical shape of the bands holds. The effect of the 
graphene lattice has been included in a recent calculation 
of the self-energy within the renormalized-ring-diagram 
approximation^. Both doped and undoped graphene 
present an imaginary part of the self-energy that, near 
the chemical potential, varies as quadratic in the energy 
concluding that electrons in graphene behave like a mod- 
erately correlated Fermi liquid^. 

The role of electron-electron interactions in ARPES 
spectra of graphene has been investigated within the 
RPA and it is shown to be important on the vkp scale^^ 
suggesting that, varying the carrier density, the effects 
of electron-electron interaction can be separated from 
electron-phonon interaction effects. The anisotropy of 
graphene energy-constant maps in ARPES has been pro- 
posed to characterize quasiparticle properties^ from the 
electronic chirality on monolayer graphene to the mag- 
nitude and sign of the interlayer coupling on bilayer 
graphene. Information about substrate-induced asymme- 
try may also be extracted from ARPES constant energy 
patterns^. 

Due to the high potential of photoemission to char- 
acterize graphene samples, the knowledge of the many- 
body effects is crucial to interpret some feature in the 



experimental spectra. The topology of the Fermi surface 
is related with fundamental properties of the material. 
To keep track of the different factors involved in the ex- 
perimental data is needed when comparing theory and 
experiment. As far as we know, corrections to the FS 
topology of 2D doped graphene, the main focus of the 
present work, have not been specifically addressed so far. 
We conclude that our results highlight the importance of 
electronic correlation in graphene. 
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APPENDIX A: COULOMB INTERACTION 
BETWEEN NEAREST NEIGHBOR ATOMS IN 
THE HONEYCOMB LATTICE. 

The Coulomb interaction between electrons in nearest- 
neighbor atoms in graphene has the form 



V =sl if 

a 





otherwise 



(Al) 



where ^ stand for lattice vectors in the two distinct 
A and B sites of the honeycomb lattice. Each atom of 
the sublattice A have three nearest neighbor atoms in the 
sublattice B, and viccversa, fj = where fj are 

vectors connecting nearest-neighbors sites on the sublat- 
tice B from the sublattice A (see Fig 1(a) I. The nearest- 



neighbor interaction term of the Hamiltonian is given by 
(we omit spin indices) 



R J = l 

where R, = riiSi -I- n2a2 is the position of an A atom 
in the lattice, with ni.2 integer numbers. The above op- 
erators in the momentum space are expressed as 



.fTT E 



^-ik-(R+f,)rt 



hi 



(A3) 



where N is the number of unit cells. The nearest- 
neighbor Coulomb interaction is found to be 
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V 



V 



(A4) 



where /(q) = X] j=i e'^ '^^ and the sublatticc density 



operators are defined as = X^kpsz -^k equiv- 



alent definition for p?). The operators given in Eq. 



are related to the creation operators in the subband basis 
by the rotation ct = l/\/2 (e^'-'^k/'^ai + Xe^'^k/'^bi] . In 
this basis, the interacting Hamiltonian can be written as 



Til 



E 

A,A'=± q£BZ 



with the expression given in Eq. (fTO|) for F*''(q). Note 
that, in addition, the overlap matrix elements of the elec- 
tronic wave-functions must be taken into account in the 
computation of the exchange self-energy. 



APPENDIX B: INTEGRAL IN MOMENTUM 
SPACE. 

In order to calculate the exchange self-energy from 
Ea. (jl3p we carry out a momentum integral in the first 



Brillouin zone. We divide the BZ in two different regions 
as shown in Fig 6(a) The long-range interaction dom- 



inates in the central region while the nearest-neighbor 
interaction is important in the outer region of the BZ. 
The circle centered in F, < |q| < 5-^^ is the boundary 
between the two regions. 

The interpolated potential considered in the calcula- 
tion is schematically plotted in Fig |6(b)[ The total po- 
tential, resulting from the combination of the long and 





(a) 



(b) 



FIG. 6: (a) Integration regions in the hexagonal Brilluoin 
zone, (b) Schematic plot of the Coulomb potential used in 
the calculation. The black circle separates the inner region 
where the long-range interaction is dominant from the outer 
region where the short range-interaction dominates. 

short range interactions, is a continuous function consis- 
tent with the lattice symmetry. 
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